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One metric for comparing confinement properties of different magnetic fusion energy configurations is tlie 
linear critical gradient of drift wave modes. The critical gradient scale length determines the ratio of the core 
to pedestal temperature when a plasma is limited to marginal stability in the plasma core. The gyrokinetic 
turbulence code GS2 was used to calculate critical temperature gradients for the linear, collisionless ion tem- 
perature gradient (ITG) mode in the National Compact Stellarator Experiment (NCSX) and a prototypical 
shaped tokamak, based on the profiles of a JET H-mode shot and the stronger shaping of ARIES- AT. While a 
concern was that the narrow cross section of NCSX at some toroidal locations would result in steep gradients 
that drive instabilities more easily, it is found that other stabilizing effects of the stellarator configuration 
offset this so that the normalized critical gradients for NCSX are competitive with or even better than for 
the tokamak. For the adiabatic ITG mode, NCSX and the tokamak had similar critical gradients, though 
beyond marginal stability, NCSX had larger growth rates. However, for the kinetic ITG mode, NCSX had a 
higher critical gradient and lower growth rates until o/Lt « 1.5 a/ Lj-^crit, when it surpassed the tokamak's. 
A discussion of the results presented with respect to a/Lr vs R/Lt is included. 



I. INTRODUCTION 



II. NCSX VS. A SHAPED TOKAMAK 
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Two of the main magnetic fusion energy designs are the 
axisymmetric tokamak and the non-axisymmetric stel- 
larator. Tokamaks have seen significant heat loss due 
to turbulenceji while stellarator losses have tradition- 
ally been dominated by their larger neoclassical trans- 
port. Studying turbulent transport in stellarators is in- 
creasingly important, however, as modern stellarator de- 
signs (such as Wendelstein 7- AS (W7-AS))^ Wendelstein 
7-X {W7-X),^-^ the National Compact Stellarator Ex- 
periment (NCSX) ,5 the Large Helical Device (LHD)S, 
and the Helically Symmetric Experiment (HSX)2r— ) have 
shown or are designed to have improved neoclassical 
confinement and stability properties. Therefore, turbu- 
lence could be increasingly relevant in stellarator exper- 
iments. Several gyrokinetic studies of drift-wave-driven 
turbulence in stellarator geometry have been done.— 
Progress has even been made in optimizing stellarator 
designs to have reduced turbulent transport 



Besides comparing good stellarator configurations, one 
would like to compare stellarators and tokamaks. These 
considerations are important when designing the next 
generation of experiments. A few previous comparison 
studies have been done, such as those in Refs. |l8j and 
[iqI . Here, the gyrokinetic turbulence code GS32S is used 
to compare microinstability of the electrostatic adiabatic 
ion temperature gradient (ITG) mode and the electro- 
static collisionless kinetic ITG in the quasi-axisymmetric 
National Compact Stellarator Experiment (NCSX) de- 
sign to that of a highly-elongated tokamak. Because their 
geometry differs so significantly, it is hard to pinpoint 
what parameter has the greatest effect, but overall ef- 
fects will be examined. 



To understand the trade-offs between stellarator and 
axisymmetric geometry and their confinement capabil- 
ities, designs can be compared computationally. One 
metric of confinement quality is the ratio of the core 
temperature to the pedestal temperature, To/Tped, as 
fusion reactors need very high core temperatures, and 
high core temperature implies good confinement. This 
ratio is related to the critical temperature gradients. If 
—dT/dr sa T / Lt.ctU, temperature-gradient-driven insta- 
bilities are marginally stable-a reasonable assumption in 
a reactor plasma, as temperatures inside the pedestal will 
be so high that profile stiffness will ensure that gradients 
are close to marginal stability. This is demonstrated by 
Fig. 3 of Ref. [2lj, which shows that fusion power (and 
thus the temperature profile) depends primarily on the 
pedestal temperature and not the beam power, for the 
case of balanced beams. 

At marginal stability (assuming l/Lr.crit is indepen- 
dent of minor radius). 



T(r) = Toe- 



-r/Lj 



(1) 



The minimum temperature is at the edge where r is max- 
imum, r„iax = o,, where a is the minor radius. In this sim- 
plified story, another approximation will be made, that 
Tped occurs at r = a. So, T{a) = Tped = Toe""/-^^ "". 
Therefore, the core temperature's dependence on the crit- 
ical temperature gradient for typical tokamak valnes^^— 
of a/i? w 1/3.5 and R/Lx.crit « 5 is 



(2) 



To/Tp,d = e"/^^.-" 

~ g(l/3.5)5 ^ 4 2 

One wants to maximize the core temperature, Tq. Tped 
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is set by non-transport mechanisms and cannot be arbi- 
trarily high, leaving a/ Lx^crit as the important parame- 
ter in equation O If an alternative fusion device design 
could increase a/LT,crit by just 30%, this would increase 
the central temperature by 50%, and more than double 
the fusion power. (A caveat-this is a simple estimate, 
and does not take into account the fact that MHD sta- 
bility changes with higher pressure peakedness.) In this 
paper, the critical temperature gradients are compared 
for NCSX and a strongly-shaped tokamak design. The 
stated stellarator minor and major radii are the average 
values. 



A. Miller equilibrium for this tokamak 
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NCSX runs were compared to a potential high- 
elongation tokamak based on a composite of ARIES- 
AT^ and JET H-mode shot #52979. It is well known 
that tokamak performance improves at high elongation 
and triangularity^^"— (in large part because this leads 
to more plasma current at fixed q), so when designing 
future tokamaks, one would like to use the highest possi- 
ble values of elongation and triangularity, though elonga- 
tion is limited by vertical stability control if it becomes 
too large. Some initial studies of shaping effects with 
GS2 were carried out in Ref. 24, using a range of shapes 
scaled from the particular JET shot #52979 (available in 
the ITER profile database^S). This JET shot, described 
in more detail in Refs. Is^ and [HI was chosen as a repre- 
sentative H-mode plasma that has been studied in detail 
by gyrokinetics codes before. Parameters for this paper 
were chosen from shot #52979, but the shaping param- 
eters were scaled to the higher levels achievable in toka- 
maks. These values of edge elongation and triangularity 
are from ARIES- AT, which generally tried to maximize 
these parameters subject to engineering and vertical sta- 
bility constraints. This JET shot has a conventional q 
profile, while the ARIES- AT design study assumed that 
a reversed shear scenario can be stably maintained in 
steady state. The composite tokamak of this paper has a 
conventional q profile. The Miller equilibrium'^^ for this 
prototypical or generic strongly-shaped tokamak was set 
up in the following way. 

The shaping study in Ref. [13 chose to focus on the 
radius r/a ~ 0.8 in order to be fairly near the plasma 
edge where shaping effects are stronger, but not too far 
out because because gyrokinetic codes often do not com- 
pare well with experiments near the edge (perhaps be- 
cause edge turbulence is driven by mechanisms other 
than ITG/TEM modes that require higher resolution 
than usual or additional effects that are not included in 
present gyrokinetic codes). Therefore, the study in this 
section uses r/a = 0.8. 

Just inside the separatrix at the 95% poloidal flux sur- 
face, the JET shot had an elongation of Kgs = 1.73 and 
triangularity of (^95 = 0.46, while the core elongation 
is Kcore ~ 1-3. At the radius of interest, r/a = 0.8, 



FIG. 1. Illustration of flux surface shapes for a prototype 
strongly-shaped tokamak at r/a = 0.8 (blue solid line), 0.9 
(green dashed line), and 0.98 (red dash-dot line), (color on- 
line) 
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jQ g — 0.60, and the Shafranov shift is dRo/dr 



= 0.57, So.8 = 0.19, 
-0.14. 

2.03 and magnetic shear 



R/a = 3.42, Kns = 1-46, 

Si 

Finally, the safety factor (70.8 
s = 1.62. 

Keeping JET's Shafranov shift, q, and s, a modified 
ARIES-AT case was created using its K95 — 2.08 and 
S95 — 0.76. Assuming, from the tokamak shaping studies, 

that k' CX (K95 — Kcore) ^ud S,S' OC 6gc^: 
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Representative flux surfaces for this prototype 
strongly-shaped tokamak are shown in Figure [TJ 



B. NCSX geometry 

Refs. [T^ and [s^ describe how non-axisymmetric ge- 
ometry input is created for GS2. The coordinate system 
of the fiux-tube code GS2 includes the radial coordinate, 
p — -y/s (s ~ (r/a)'^ is the normalized toroidal fiux), the 
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Parameter 


Value 


r /a 


0.8 


s^{{r/a)f 


0.64 


a = C~qe 





eo 







1.70 
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0.835 




0.0% 


R 


~ 4.7ajv ~ 1.5m 


an 


^ 0.322m 


Ba = (B) 


i.58r 



TABLE I. Geometry values for the NCSX equilibrium. 



coordinate aligned to the field line, 9, and the angle that 
selects a flux tube, a = ( ~ q{9 — 9o) (where C and 9 are 
Boozer toroidal and Boozer poloidal coordinates and 6*0 
is the ballooning parameter). 

The GS2 documentation'^^ defines geometrical quanti- 
ties in terms of a parameter d^* n / dp, where p is the radial 
coordinate and ^Pat is the normalized poloidal flux. Geo- 
metrical quantities in this paper follow GS2 notation and 
include d'^N /dp. For more information, see Refs. ITBI and 

Figures (2111] show the magnitude of the magnetic field, 
curvature drift, and for both the tokamak and 

NCSX field lines, for the entire domain and a close-up 
around 9 = Q. See Table U for a complete list of geomet- 
rical quantities and their values. 

Notice that the bad (positive) curvature regions of 
NCSX are much more localized than the tokamak case. 
Coupled with the much stronger local magnetic shear 
(responsible for the sharp peaks in fcj_ oc s. Fig. IH]), this 
explains why NCSX's electrostatic potential eigenfunc- 
tions are also more localized than the tokamak's. An 
example is shown in Figure [51 In addition, these traits 
could predict better transport properties for NCSX. 
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FIG. 2. The NCSX (blue solid hne) and tokamak (green 
dashed line) equilibria: normalized \B\ vs. 9. (color online) 




FIG. 3. The NCSX (blue solid line) and tokamak (green 
dashed line) equilibria: normalized \B\ vs. 9, showing a close- 
up around 9 — 0. (color online) 



C. ITG mode with adiabatic electrons 

For the initial study, the ITG mode with adiabatic 
electrons growth rates and their dependence on temper- 
ature gradient were compared. Figure [MTUl show typical 
growth rate spectra for NCSX and this tokamak. In Fig- 
ure [Til the growth rate at each a/Lx {a/Lxe = o-/LTi) 
was the highest in the range kyPi e [0.2,1.4] for NCSX 
and kyPi G [0.1, 1.0] for the tokamak. These ranges were 
wide enough to capture the peak of the growth rate spec- 
trum. Growth rates shown are normalized such that 
(7,^) = {■Jphysicai,ujphysicai){a/vthi)- The NCSX thresh- 
old is a/LT,crit ~ 1-26 and the tokamak's is a/Lx^crit ~ 
1.22. This difference is not very significant. However, 
soon after the threshold, the NCSX growth rates surpass 
those of the tokamak, indicating that for a given a/Lx, 
the adiabatic ITG mode is more unstable in NCSX than 



in the tokamak. This implies that the transport would 
be stiffer, but the temperature gradients would still be 
expected to be very similar since they would be set by 



D. ITG mode with kinetic electrons 

The threshold of the ITG mode with kinetic electrons 
(with a/Ln = 0) for the tokamak was somewhat lower 
than that of NCSX, but the slope of the growth-rate 
curve is almost the same for both (Fig. [T^ . Similar 
to Section III C[ growth rates shown were the highest on 
a spectrum of kyPi e [0.2, 1.4] for NCSX and kyPi e 
[0.1, 1.0] for the tokamak (see Figure [T51im show typical 
growth rate spectra). With kinetic electrons, the growth 
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FIG. 4. The NCSX (blue solid line) and tokamak (green 
dashed line) the curvature drift frequency {iOcv,norm ~ 
(2a^/Bjv)(d*jv/dp)(fcx/n)-bx[b-Vb]) along 6* . (color online) 
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FIG. 6. The NCSX (blue solid hue) and tokamak (green 
dashed line) equilibria: ^■^^ '"s. 9. (color online) 




FIG. 5. The NCSX (blue sohd hue) and tokamak 
(green dashed line) equilibria: the curvature drift frequency 
(i^cv.norm = {2a^ / B n) (d'if N / dp)(k ± / u) ■ b X [b ■ Vb]) along 
9, showing a close-up around 9 = 0. (color online) 




FIG. 7. The NCSX (blue solid hue) and tokamak (green 
dashed line) equilibria: {^1^^ ^' showing a close-up 

around 9 = 0. (color online) 



rates for the ITG mode in NCSX increased over the adi- 
abatic electron case (Fig. [TT|) . while the critical gradi- 
ent lowered to a/LT,crit ~ 1.21. The tokamak threshold 
decreased somewhat further, to a/LT,crit ~ 1.11. The 
slope of the NCSX line is somewhat steeper, and for 
q/Lt ~ 1.82, the NCSX growth rates are larger than 
the tokamak growth rates. 

The growth rate vs. a/Lx plot in Figure [T^ shows im- 
provement in the a/Lx threshold of NCSX over this toka- 
mak by about 10%. Based on the marginal stability logic 
in the beginning of Section [III this corresponds to about 
22% more fusion power for a NCSX-based design rela- 
tive to the tokamak (with the same edge temperature 
and density assumed for the two designs, and approxi- 



mating the fusion power as scaling as T^). (Other effects 
that might change this result include finite beta modifi- 
cations to the equilibrium and, nonlinearly, the Dimits 

Shift.22^) 

This is much better than one might have initially 
guessed based on just the local value of R/Lt in NCSX 
vs. a tokamak. While from equation [5] it is clear that 
a/ Lt is the relevant parameter for determining the core 
temperature, in the axisymmetric community, the thresh- 
old for the ITG instability is usually expressed in terms 
of R/Lt^ which is often the key parameter numerically. 
An instability threshold R/LT.crit can be derived from 
the dispersion relation for a local ITG mode in the bad 
curvature region, ignoring the parallel dynamics. In 
this limit, the critical instability parameter is the ra- 
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FIG. 8. Comparing electrostatic eigenfunctions for NCSX 
{Re{(j}): red triangles and Im{4>): light blue solid line) and 
tokamak {Re{())): blue circles and Im{(f)): green dashed line), 
for an adiabatic ITG mode with a/Lr = 'i,a/Ln — 0. For 
ARIES, kypi = 0.55, and for NCSX, kypt = 1.0. (color online) 
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FIG. 9. NCSX growth rate spectrum for the adiabatic ITG 
mode with o/Lt ~ 3,a/L„ — 0. (color online) 
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FIG. 10. Tokamak growth rate spectrum for the adiabatic 
ITG mode with a/Lr = 3, a/L„ — 0. (color online) 
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FIG. 11. NCSX (blue crosses) and ARIES- AT-like tokamak 
(red circles) adiabatic ITG mode growth rate dependence on 
temperature gradient. Fits obtained through piecewise lin- 
ear interpolation on the lowest half of the growth rate curve, 
(color online) 



tio of the temperature-gradient diamagnetic drift fre- 
quency (a;*T oc l/Lxi) to the curvature drift frequency 
{ujd (X 1/i?) . (Particles with different energies have differ- 
ent curvature drift velocities, which can result in Landau 
damping. This criterion essentially says that the drive 
from the temperature gradient must be strong enough to 
overcome this damping in order to drive instabilities.) 

A concern could be that if an NCSX design is limited 
to the same local Rioc/ Lt,Ioc sls in a tokamak, it would 
have a much lower g/Lt (because at some toroidal lo- 
cations, such as the left panel of Figure [151 the cross 
section of NCSX is very narrow, with a local plasma half- 
width aioc ~ a/2.58), and thus would have much lower 
fusion power. Therefore, the local value of the logarith- 



mic gradient, 1/Lt.Ioc = |Vr|/r = {a/ LT)/aioc is much 
larger than the average 1/Lt- This is enhanced by the 
larger average aspect ratio {R/ajMCSX = 4.7 relative to 
the tokamak {R/a)tok — 3.42, and is partially compen- 
sated by the fact that the local radius of curvature of 
the magnetic field, i?/oc = \i> ■ V6|^^ = 0.92m (evaluated 
at the outer midplane of the r/a = 0.8 fltix surface in 
the left panel of Figure [TS]) is somewhat smaller than 
the average radius of curvature R = 1.51m in NCSX. 
Considering these modifications, the local Rioc/ Lt loc = 
{Rioc/R){alaioc){R/a){a/LT) = 7.39 a/Lr in NCSX, 
while R/Lt = 3.42 a/Lx for a tokamak. 

Restating this concern, one may have thought that if 
NCSX and a tokamak had the same normalized temper- 
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FIG. 12. Growth rates for an ITG mode with kinetic electrons 
as a function of temperature gradient for NCSX (blue crosses) 
and an ARIES- AT- like tokamak configuration (red circles). 
Fits obtained through piecewise linear interpolation on the 
lowest half of the growth rate curve, (color online) 
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FIG. 13. NCSX growth rate spectrum for the kinetic ITG 
mode with a/Lr ~ 3, a/L„ — 0. (color online) 
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FIG. 14. Tokamak growth rate spectrum for the kinetic ITG 
mode with a/Lr = 3, a/Ln = 0. (color online) 
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FIG. 15. Poloidal cross sections of NCSX for two toroidal 
angles. The dashed line is the location of the vacuum vessel 
and the solid lines are last closed flux surfaces for various t 
profiles. Figure 2 of Ref. [35l : reprinted with permission. More 
information can be found in Ref. [3q. 



bad-curvature region (i.e., a shorter connection length 
between good and bad curvature regions) and stronger 
local magnetic shear. 



ature gradient, a/Lx, the ITG modes would be much 
worse in NCSX due to a much higher Rioc/ Lt,Ioc than 
the tokamak. In fact. Figure [T2l showed that NCSX has 
a somewhat higher critical gradient in terms of a/Lx, 
so the hypothesis that NCSX and a tokamak are similar 
when expressed in terms of Rioc/ Lt,Ioc must be incorrect. 
Indeed, this is strikingly illustrated in Figure [16] (same 
data as Fig. [121 renormalized) , which shows that NCSX 
has in fact much lower growth rates than a tokamak for 
the same Rioc/ Lt,Ioc- This is probably because the par- 
allel dynamics are in fact not negligible in NCSX. The 
eigenfunctions, as seen in Figure [51 are more localized 
along a field line in NCSX than in a tokamak, through 
some combination of the stabilizing effects of a narrower 



III. CONCLUSION 

Cross-configuration comparisons of plasma confine- 
ment are important to consider for the design of future 
fusion energy devices. As a simple case, the linear stabil- 
ity was compared for NCSX and a tokamak equilibrium. 
This particular tokamak equilibrium is a composite of 
JET H-mode shot #52979 and ARIES- AT. NCSX had a 
similar linear critical temperature gradient a/Lx^crit to 
the tokamak case for ITG modes with adiabatic electrons, 
though its growth rates were higher than the tokamak's 
beyond marginal stability. However, for ITG modes with 
kinetic electrons, NCSX's critical gradient u/Lt is ap- 
proximately 9% higher than the tokamak's, which would 
correspond to an approximately 22% increase in the fu- 
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